The Classical-Map Hyper-Netted-Chain (CHNC) method and associated novel 
density-functional methods for Warm Dense Matter* 
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he advent of short-pulse lasers, nanotechnology, as well as shock-wave techniques have created 
new states of matter (e.g., warm dense matter) that call for new theoretical tools. Ion correlations, 
electron correlations as well as bound states, continuum states, partial degeneracies and quasi- 
equilibrium systems need to be addressed. Bogoliubov's ideas of timescales can be used to discuss 
the quasi-thermodynamics of non-equilibrium systems. A rigorous approach to the associated many- 
body problem turns out to be the computation of the underlying pair-distribution functions gee, Qei 
and ga, that directly yield non-local exchange-correlation potentials, free energies etc., valid within 
the timescales of each evolving system. An accurate classical map of the strongly-quantum uniform 
electron-gas problem given by Dharma-wardana and Perrot is reviewed. This replaces the quantum 
electrons at T = by an equivalent classical fluid at a flnite temperature Tg, and having the same 
correlation energy. It has been shown, but not proven, that the classical fluid Qij are excellent 
approximations to the quantum Qij. The classical map is used with classical molecular dynam- 
ics (CMMD) or hyper-netted-chain integral equations (CHNC) to determine the pair-distribution 
functions (PDFs), and hence their thermodynamic and linear transport properties. The CHNC is 
very eflicient for calculating the PDFs of uniform systems, while CMMD is more adapted to non- 
uniform systems. Applications to 2D and 3D quantum fluids. Si metal-oxide-field-effect transistors, 
Al plasmas, shock-compressed deuterium, two-temperature plasmas, pseudopotentials, as well as 
calculations for parabolic quantum dots are reviewed. 

PACS numbers: PACS Numbers: 71.10.Lp,75.70.Ak,73.22-f 



MOTIVATION 

The advent of powerful short-pulse lasers as well 
as other new tools for manipulating matter presents 
new challenges to existing theory. Warm dense matter 
(WDM) is such a regime where we have highly corre- 
lated ions, electrons, finite temperature as well as partial 
degeneracy effects that have to be taken into account. 
Sufficiently thin nano-slabs of WDM can be studies with 
a variety of probes [H 0] . Bound states as well as contin- 
uum states have to be treated without sinking in a morass 
of computations. The Born-Oppenheimer approximation 
cannot be used if coupled-mode effects are important. In 
this paper we examine new theoretical approaches that 
extend beyond the familiar territory of density- functional 
theory (DFT) to treat these and other intractable prob- 
lems in many-body physics. 

The Hohenberg-Kohn and Mermin (HKM) theo- 
rems Q of DFT assert that the one-body density n{r) 
of an inhomogeneous system completely determines its 
physics. However, implementations of DFT use the more 
laborious Kohn-Sham (K-S) approach Q in lieu of an ac- 
curate kinetic-energy functional The Kohn-Sham 
n{r) of an electron system is: 



The K-S eigenstates, ipi, with "energies" e^, occupations 
factors /i/ at the temperature T = for all the quan- 
tum numbers v have to be determined, self-consistently, 
using a one-body Kohn-Sham potential Vks in the Kohn- 
Sham equation. The inclusion of continuum states in 
this summation consistently, to satisfy sum rules etc., is 
a challenge discussed in [7|. The Kohn-Sham potential 
contains an 'exchange-correlation potential' 14;c(['t-]) that 
maps the many-body effects to a functional of the one- 
body density. Model Fa;c([n]) potentials have been con- 
structed using microscopic theories of systems like the 
uniform electron liquid (UEF). Such UEF-calculations 
are equivalent to a coupling-constant integration over the 
electron-electron pair distribution function (PDF), viz., 
gee(r)- Calculating these PDFS, even for uniform sys- 
tems, is a challenge that is treated in this paper. 

Quantum systems at high temperatures behave clas- 
sically. Then the Kohn-Sham procedure simplifies. The 
density n(f) is given by the Boltzmann form: 



2(f) = no exp{-^VR-s(^)} 



(2) 



(1) 
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where no is a reference density, and Vks is a classical 
Kohn-Sham potential that has to be obtained from a mi- 
croscopic classical many-body theory. The 'potential of 
mean-force' used in classical liquid-state theory is just 
this classical Vks- If the center of coordinates is selected 
to be one of the classical particles, and if we consider a 
uniform fiuid, then n{r) becomes the density profile of 
field particles around the central particle which acts like 
an external potential. The density profile n(r) is directly 
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related to the pair-distribution function, i.e., 
n{r) — ng{r), n — n{r — > oo). 



(3) 



Hence one may attempt to go beyond traditional DFT 
and proceed directly to the underlying calculation of 
the pair-densities themselves. The extension of the 
Hohenberg-Kohn theorem given by Gilbert, using the 
one-body reduced density matrix is actually entirely in 
this spirit Q. However, the PDF is conceptually eas- 
ier to use than the density matrix. Such considerations 
suggest that the kinetic-energy functional may be side- 
stepped by: (i) the use of an equivalent "classical-fluid" 
at a temperature Tcf for the (uniform) quantum fluid 
whose actual physical temperature T may even be zero; 
(ii) the use of effective classical pair-potentials inclusive 
of quantum effects to calculate classical pair distribution 
functions which can then be used to compute most of the 
usual physical properties 13[. 

The advantage of such a classical-map approach is 
that the ions, being essentially classical particles, can be 
treated together with the electrons in the same classi- 
cal computational scheme. Unlike quantum iV-electron 
schemes which, in principle, grow in complexity non- 
polynomially in iV, classical methods are essentially inde- 
pendent of N . Here we should note that traditional quan- 
tum chemistry and condensed-matter physics treat only 
the electrons by DFT. On the other hand. Gross and col- 
laborators have attempted to present a completely quan- 
tum mechanical non-adiabatic DFT theory of electron- 
nuclear systems, and given an application to the H2 sys- 
tem P|. In standard calculations, the ion positions are 
explicitly included and form the external potential for 
the motion of the Kohn-Sham electron. In warm dense 
matter (WDM), e.g., highly compressed hot hydrogen, 
there are as many protons as there are electrons in a 
given volume of the sample. Ion motion couples with 
electron-plasma oscillations to generate ion-acoustic cou- 
pled modes. Their effects may be missed out in standard 
DFT formulations as well as in MD simulations. 

In any case, the quantum-chemistry approach (e.g., as 
in the Gaussian package) rapidly becomes intractable, 
esp. when continuum states have to be included - as in a 
plasma. The solid-state approach of using a periodic cell 
is more flexible here, as in the Vienna-simulation pack- 
age known as the VASP. However, WDM applications de- 
mand large unit cells and calculations of energy bands for 
many ionic configurations. The classical-map approach, 
where both ions and electrons are treated as classical flu- 
ids inclusive of particle motions, provides a new paradigm 
for warm dense matter and other novel systems which are 
computationally very demanding by standard methods. 
Such standard methods could be regarded as microscopic 
bench marks for more global methods like the classical- 
map approach discussed here. 

The philosophy of the classical-map technique is to 
treat the zeroth-order Hamiltonian exactly, i.e., using 



the known quantum solution, and then use the classical 
map for dealing with the many-body effects generated 
from the Coulomb interaction. For uniform systems, the 
eigen-solutions of the Hq problem are plane waves. Fermi 
statistics imposes a determinantal form to the wavefunc- 
tions, and hence the non-interacting PDFs g^^, are differ- 
ent from unity if the spin indices s, s' are identical. Thus 
g^s{r) exhibits a Fermi hole, which can be exactly rep- 
resented by a classical repulsive potential known as the 
Pauli exclusion potential (PEP) This should per- 

haps be called the 'Fermi-hole potential' as it should not 
be confused with the 'Pauli Potential' defined in DFT 



111 |12| via the density-functional derivative of the dif- 
ference between the non-interacting kinetic energy and 
the full von Weizsacker kinetic energy. In the interest of 
historical accuracy, it should however be noted that the 
name 'Pauli potential' was already in use for the Fermi- 
hole potential since the work of Lado. We use the names 
'Fermi-hole potential' and Pauli-exclusion potential' as 
synonymous, and different from the DFT correction to 
the von Weizsacker term known as the Pauli potential. 
The key ingredients of the method are the following. 

1 . Replacement of the electron system at T by a classi- 
cal Coulomb fluid at an effective classical-fluid tem- 
perature given by 



T,f = (T^ + T 



n2\l/2 



(4) 



where Tq is a 'quantum temperature' which de- 
pends only on the electron density. Tq is such that 
the classical fluid has the same correlation energy 
as the initial quantum fluid at T = 0. The motiva- 
tion for defining Tcf by Eq. H] is given in [iJl • 

2. Inclusion of a "Pauli exclusion potential", i.e., a 
Fermi- hole potential (FHP) to reproduce the Fermi 
hole of spin-parallel electrons exactly. 

3. The use of a diffraction-corrected Coulomb interac- 
tion (l/r)(l — e"''/'*') to account for the finite-size 
of the de Broglie thermal wavelength A of the elec- 
trons at the finite temperature Tq. 

4. Calculation of the pair-distribution functions of the 
classical fluid using an integral-equation method 
(CHNC), or molecular dynamics. When MD is 
used in this manner we call it classical-map molec- 
ular dynamics (CMMD). 

5. use of the PDFs in coupling-constant integrations 
to calculate the Helmholtz free energy and all other 
thermodynamic properties of the quantum fluid. 
The linear transport properties (e.g., conductivity) 
are available from Kubo or Ziman-type formula- 
tions which use the PDFs and potentials as inputs. 

Formulations which use this method have been success- 
fully applied to a number of quantum systems: 
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(i) The 3-D electron fluid at T = and at finite T [13 1 



(ii) The 2-D electron fluid both at T = [ISHlSl, and at 
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finite T 

(iii) The calculation of Fermi-liquid properties like the 
electron effective mass m* , the enhancement of the Lande 
(/-factor IiqII . and local-field corrections to the response 
functions '20|. 

(iv) The multi-component electron fluid in Si-Si02 metal- 
oxide- field-effect transistors [2l|; preliminary applica- 
tiOTLS to multi-valley massless Dirac fermions in graphene 

(v) Electrons confined in parabolic potentials (quantum 
dots) 

(vi) Two-mass two-temperature plasmas [25| . 

(vii) Equation of state and Hugoniot of Shock-compress 
hydrogen [26| . 

(viil) Liquid Al under WDM conditions; linear transport 
properties of some WDM systems, where some of the 
PDFs were calculated using CHNC [27 j. 



The QHNC method of Chihara 

For the sake of completeness we also mention Chi- 
hara's 'quantal-HNC (QHNC) method @. Here an 
HNC-type equation is solved for the electron subsys- 
tem. The electron-electron pair-distribution function is 
calculated by solving the 'quantal HNC equation' with 
"a fixed electron" at the origin. However, the electron 
pair-distribution functions obtained by this method for 
jellium are in poor agreement with those from quantum 
Monte Carlo methods. In fact, if non-interacting elec- 
trons are considered, the zeroth order PDF, which is 
known analytically at T = and in terms of a Fermi 
integral at finite-T (as discussed below) is not recovered 
correctly by Chihara's method. The small-fc limit of the 
ion-ion structure factors calculated by QHNC fail to re- 
produce the correct compressibility. Nevertheless, Chi- 
hara's QHNC recovers some of the short-ranged order 
in the ion-ion pair-distribution functions, where the os- 
cillations and peak heights are in rough agreement with 
microscopic simulations. The short-comings in Chira's 
formulation are overcome in the CHNC method. 

In the following we discuss details of some of the im- 
plementations of CHNC using integral-equation meth- 
ods since they are conceptually more transparent and 
far cheaper than molecular dynamics (CMMD), let alone 
QMC. 



A CLASSICAL REPRESENTATION FOR THE 
UNIFORM ELECTRON LIQUID 

A system of electrons held in place by an external po- 
tential (as in a solid, a quantum well, or in a molecule) 
at T = is necessarily a quantum system. The uniform 



electron fluid (UEF) at a density n, Wigner-Seitz radius 
r^, is the key paradigm for treating exchange and corre- 
lation in DFT. The pair-distribution functions (PDFs) of 
the UEF at T = are known from quantum-Monte Carlo 
(QMC) studies. They are the basis of exchange-and cor- 
relation energies of the UEF. Hence, if the classical-map 
scheme could successfully calculate the PDFs of the elec- 
tron fluid at arbitrary coupling and spin polarization, in 
2-D and 3-D, then the idea that the quantum fluid can be 
represented by a classical Coulomb fluid stands justified. 

Consider a fluid of mean density n containing two spin 
species with concentrations Xi = rii/n. We deal with the 
physical temperature T of the UEF, while the tempera- 
ture Tcf of the classical fluid is 1//3. Since the leading 
dependence of the energy on temperature is quadratic, 
we construct Tcf as in Eq. This is clearly valid for 
T = and for high T. This assumption has been exam- 
ined in greater detail by various applications where it has 
been found successful. 

The properties of classical fluids interacting via pair 
potentials (pij (r) can be calculated using classical molec- 
ular dynamics (MD) or using an integral equation like 
the modifled hyper-netted-chain equation. The pair- 
distribution functions for a classical fluid at an inverse 
temperature /3 can be written as 

g,j (r) = exp[-l3(l),j (r) + h,j (r) - Cy (r) + Bij (r)] (5) 

Here 4>ij{r) is the pair potential between the species 
For two electrons this is just the Coulomb poten- 
tial Vcoui'"')- If the spins are parallel, the Pauli exclusion 
principle prevents them from occupying the same spa- 
tial orbital. Following the earlier work, notably by Lado 
[lot , we also introduce a "Pauli exclusion potential" or 
Fermi- hole potential (FHP), V{r). Thus (t>ij{r) becomes 
V{r)Sij + Vcou{r). The FHP, r{r), is constructed to re- 
cover the PDFs of the non-interacting UEF, i.e., gij{r) is 
exactly recovered. The function h(r) = g{r) - 1; it is re- 
lated to the structure factor S{k) by a Fourier transform. 
The c(r) is the "direct correlation function (DCF)" of the 
Ornstein-Zernike (OZ) equations. 



%(r) 



(|r-r'|)c,,,(r') (6) 



The Bij (r) term in Eq.[5]is the "bridge" term arising from 
certain cluster interactions. If this is neglected Eqs. [Sjn] 
form a closed set providing the HNC approximation to 
the PDF of a classical fluid. Since the cluster terms be- 
yond the HNC approximation are difficult to calculate, 
they have been modeled approximately using the theory 
of hard-sphere liquids (29i]. We have provided explicit 
B{r) functions for the 2-D electron fluid where it is im- 
portant even at low coupling (soj . B{r) is important in 
3-D when the coupling constant F = P/ts for electron- 
electron interactions exceeds, say, 20. The range of F 
relevant to most WDM work (e.g., F ~ 4.5 even for = 
10 ) is such that the HNC-approximation holds well. 
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FIG. 1; (Online color) The Ferm-hole potential (i.e., Pauli- 
exclusion potential) P'P{r) is a universal function of rkp at 
each T and reproduces the Fermi hole in the parallel-spin 
zeroth-order PDF, gss{i~), shown in the inset, for T = and 
T/Ef ~ 2. If the spins are anti-parallel, flags' (?■) = 1 and the 
Pauli-exclusion potential is zero. 

Consider the non-interacting system at temperature T, 
with Xi = 0.5 for the paramagnetic case. The parallel- 
spin PDF, i.e, 3° (r, T), wiU be denoted by (7^(r) for sim- 
plicity, since 3° (r, T), i ^ j is unity. Denoting (ri — r2) 
by r, it is easy to show, as in sec. 5.1 of Mahan [3]J, that: 

5° (r) = -lsk„k.n(A:i)n(fc2)[l - e'^''-''^)--]. (7) 

Here n{k) is the Fermi occupation number at the tem- 
perature T. Eq. [7] reduces to: 

gUr) = l-F^r) (8) 
Fr{r) = (67rV4) / n{kf-^^. (9) 

Here kp is the Fermi momentum. Thus gj^ir) is obtained 
from the Fourier transform of the Fermi function. The 
zeroth-order PDF is a universal function of rkp- It is 
shown in the inset to Fig. [TJ Assuming that g^^ (r) can 
be modeled by an HNC fluid with the pair interaction 
PVir), the "Fermi-hole potential", viz., 7^(r), is easily 
seen to be given by 

[3V{r) = -log[g"{r)] + h"{r) - c°(r) (10) 



The c°(r) can be evaluated from g^{r) using the OZ re- 
lations. The T = case can be evaluated analytically 

We can determine only the product l3V{r). The classi- 
cal fluid "temperature" 1 //3 is still undefined and clearly 
cannot be the thermodynamic temperature T as T — >■ 0. 
The Pauli-exclusion potential, i.e., the FHP, is a universal 
function of rkp at each T. It is long ranged and mim- 
ics the exclusion effects of Fermi statistics that produces 
quantum entanglement. At finite T the range of the 
Pauli-exclusion potential is comparable to the de Broglie 
thermal wavelength and is increasingly hard-sphere like. 
Plots of /3V{r) and gssi^) are given in Fig. [1] 

The next step in the CHNC method is to use the full 
pair-potential (j^ijir), and solve the coupled HNC and OZ 
equations for the binary (up, and down spins) interacting 
fiuid. For the paramagnetic case. Hi — n/2, we have: 

h,j{q) = ^K,{r) (12) 
hii{q) = cii{q) + {n/2)[cii{q)hii{q) + Ci2{q)h2iiq)] 
hi2{q) = Ci2{q) + {n/2)[cii{q)hi2{q) + Ci2(g')/i22(g)] 

(13) 

The Coulomb potential Vcoui^) needs some discussion. 
For two point-charge electrons this is l/r. However, de- 
pending on the temperature T, an electron is localized 
to within a thermal wavelength. Thus, following ear- 
lier work, e.g., Morita, and Minoo et al. 32], we use a 
"diffraction-corrected" form: 

Vcou{r) = (l/r)[l-e-'-/^-]; A*^ = {2^mT,fr^'\ (14) 

Here m is the reduced mass of the electron pair, i.e., 
m*{rs)/2 a.u., where m*(rs) is the electron effective 
mass. It is weakly dependent, e.g, ~0.96 for = 
1. In this work we take m*=l. The "diffraction correc- 
tion" ensures the correct behaviour of gi2{r —>■ 0) for all 

In solving the above equations for a given and at 
T=0, we have Tcf=Tq. A trial Tg is adjusted to obtain 
an Ec{Tq) equal to the known paramagnetic Ecirg) at 
each Tg, via a coupling constant integration. 

„ N dX f Anr'^dr , , , , , , ^ , ^ 

EUT,) = YJ -^ihn{r,X) + h,2{r,X)] (15) 

{Ex alone is obtained if A is fixed at 0). The resulting 
"quantum" temperatures Tq could be fitted to the form: 

Tq/ Ef = 1.0/ {a + b^rs + crs) (16) 

We have also presented a fit to the Tq of the 2-D electron 
system, and discussed how the 2-D and 3-D fits could be 
related by a dimensional argument. Bulutay and Tanatar 
have also examined the CHNC method, and provided fits 
to the Tq of the 2-D electron fluid [11]. 
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FIG. 2: (Online colour) (a) The g{r) from CHNC (solid lines) 
are compared with those of Dandrea et al. [13], (dashed lines) 
at rs—2. The temperature t — T/Ep. Panel (b) — 5, 
CHNC (sohd hues) for t=0 snd t=5. The g{r) at t = from 
Tanaka and Ichimaru [SJ] (dashed line), and from DMC-QMC 
[35>] (boxes), are also shown. Panel (c) rs = 5, CHNC (solid 
lines) for t=0 snd t—20. The g{r) of Tanaka and Ichimaru 
[3^ (dashed line) is also shown for t=0. 



For any given r^, given the Tq from the paramag- 
netic case, we can obtain gij{r) and £'xc(rs, C, T') [13j . 
at arbitrary unexplored values of spin-polarization C 
by solving the coupled HNC equations, or doing an 
MD calculation using the Fermi-hole potential and the 
diffraction-corrected Coulomb potential. Many analytic 
theories of electron fluids, e.g., those of Singwi, Tosi, 
Land and Sjolander i31j, Tanaka and Ichimaru, predict 
g{r) which become negative for some values of r even 
for moderate Tg. The PDFs obtained from the HNC- 
procedure are positive definite at all r^. In Fig.[2]we show 
typical results for gij (r) and comparisons with QMC- 
simulations. Our results are in excellent agreement with 
the DMC resuhs [H. 

The Tq determined from the unpolarized Ec is used 
to calculate Ec{tsX^T^) ^iiy C- The QMC results 
for £'c(''s,C) at T = agree with ours, since our gij{r) 
agree with those from MC. For example, at = 10, the 
spin-polarized —Ec is: Ceperley- Alder, 0.0209 Ry; Ortiz- 
Ballone, 0.0206 Ry fsS]; our method (CHNC), 0.0201 Ry; 




r/r 

s 

FIG. 3: (Colour online)The QMC pair-distribution function 
of a fully spin polarized (<^ = 1) 2D electron fluid at — 10, 
and T = is compared with those calculated from CHNC 
using a hard-sphere bridge function and with no bridge func- 
tion what so ever. An essentially exact fit with QMC can be 
obtained using a Coulomb bridge function [3^. 

Kallio and Piilo, 0.0171 Ry [3^. 

Most of the recent work using CHNC has been on the 
2D-electron fiuid owing to its accrued interest in nanos- 
tructures and technological applications. The electron- 
electrons interactions are stronger in reduced dimensions, 
and the use of a bridge function to supplement the CHNC 
equation is essential for accurate work (soj . However, 
even the appropriately chosen hard-disc bridge works 
quite well, as seen in Fig. [S] 

Fermi-liquid parameters of electron fluids 

It is in fact possible that in some circumstances, WDM 
may fall into the category of a Fermi liquid. Highly 
compressed electron systems have correspondingly high 
Fermi energies and hence may have a physical tempera- 
ture T « Ep. In any case, we review the calculation of 
Fermi-liquid parameters as it is an important aspect of 
the capability of a classical map to extract results in the 
strong quantum domain. 

Microscopic many-body physics allows one to calcu- 
late various quantities like the effective mass m* or the 
Lande g-factor that enter into Landau's theory of Fermi 
liquids. One would perhaps assume that a classical rep- 
resentation of a Fermi liquid would hardly be successful 
in attacking such problems. For instance, m* is usually 
calculated from the solutions of the Dyson equation for 
the one-particle interacting Green's function. If the real 
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FIG. 4: (Colour online) The effective mass m* of an ideal 
2D layer (zero thickness) obtained from CHNC are compared 
with the the Quantum Monte Carlo data of Ref. ,38.1 and 
the perturbation theory calculations of Zhang et aZ. [Syl, and 
Asgari et aZ.js^l, i.e., their calculation labeled G+&cG- / D. 
The CHNC calculation is sensitive to the choice of the bridge 
function, while the Green's-function methods are sensitive to 
the truncations used. 



part of the retarded self-energy is Si(fc,a;), the Landau 
quasi-particle excitation energy EQp{k), measured with 
respect to the chemical potential is used in calculating 
the effective mass m*. 



EQp{k) 



efe 
1 



dEQpjk )^ 
kp dk 



I k—kf 



(17) 
(18) 

(19) 



This is a very arduous calculation, and there are tech- 
nical questions about the difficulties of satisfying sum 
rules. Ward identities etc., when the Dyson equation 
is truncated in some approximation. The values of m* 
calculated by different authors using different perturba- 
tion expansions differ significantly, and from QMC results 

Mm 

By contrast, the calculation of m*, and also g* using 
CHNC is very simple because it can evaluate the free en- 
ergy F of the electron fluid as a function of the physical 
temperature T as well as the spin polarization C. The 
ratio of the interacting and non-interacting specific heats 
provides a simple evaluation of the m* , while the ratio of 



the interacting and non-interacting susceptibilities, de- 
termined from the second derivative (with respect to 
of the exchange-correlation correction to the free energy 
provides the product m*g* [l^ . 



m 



C^ICl = 1 -f 
XpIXs = 1 + 



[a2Fo(t)/9t2] 

[a^f..(C)/ac^] 

[92Fo(C)/9C'] 



(20) 
(21) 



Detailed calculations for ideal 2-D electron layers (see 
Fig. HJ, thick layers as well as for multi- valley systems 
using the CHNC method have been presented in our pub- 
lications flff|. Calculations of m* and g* for the 3-D elec- 
tron liquid using CHNC have not yet been undertaken, 
while RPA results have been given by Rice [sij . 



DENSE HYDROGEN AND OTHER PLASMAS 

Dense hydrogen, or any other fully ionized plasma is a 
direct generalization of the uniform electron-fluid prob- 
lem to include an additional component (e.g., protons), 
while removing the positive neutralizing background. Let 
us consider a fully ionized plasma with ions of charge 
Z, and density p. Then the electron density n = Zp, 
and we assume that both subsystems are at the same 
physical temperature T . The electron subsystem will 
have to be calculated at a classical-fluid temperature 
Tcf — \/{T^ + T"^) and the electron-electron interactions 
have to be diffraction corrected. On the other hand, the 
ions are classical particles and the simulations (or integral 
equations) for the ions will use the physical temperature 
T. The quantum correction Tq can be neglected for ion, 
as discussed in [2^ . The total Hamiltonian now contains 
the three terms. Hi, He, and the electron-ion interaction 
Hei- The electron system contains two spin components, 
while the ion system adds another component. Thus, a 
three-component problem involving six pair-distribution 
functions have to be calculated. If spin effects could be 
neglected, then the two spin-components of the electrons 
could be replaced by an effective one-component electron 
fluid where the Pauli-exclusion potential (i.e, FHP) is in- 
cluded after averaging over the two components. 

An example of a classical-map calculation of the EOS 
of laser-shock compressed hydrogen has been given by 
Dharma-wardana and Perrot [l^ , where a Hugoniot has 
been calculated and compared with those from other 
methods (see Fig. [5]). The article by Michael Desjarlais 
in this issue also refers to the problem of the equation 
of state (EOS) of highly compressed hydrogen [40|. A 
proper experimental probe of such laser-compression ex- 
periments needs to address some method of independent 
measurement of the electron temperature and the ion 
temperature Ti. If the electrons and ions are in equilib- 
rium, T = Ti = Tf.. Then the usual DFT methods using 
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the Born-Oppenheimer decoupling would be expected to 
give a good prediction of the EOS, and also the Hugo- 
niot. The EOS calculation is essentially a calculation of 
the partition function. This requires the evaluation of 
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< g-(^fe/T,+H,/T,+H„/T„) ^ 



(22) 



Here the total Hamiltonian H is rewritten in terms of 
iJe, and the electron- ion interaction i/gi which is 
again a Coulomb potential. We have included a cross- 
subsystem temperature T^i which is simply T for equilib- 
rium systems. If a Born-Oppenheimer approximation is 
used, the electrons 'do not know' the temperature of the 
ions, and vice versa. For equilibrium systems, a Born- 
Oppenheimer correction can be introduced, e.g., as in 
Morales et al. [41|. However, the add on correction in- 
troduced by Morales et al. will change the virial com- 
pressibility, leaving the small fc behaviour of the 
proton-proton structure factor unaffected, and hence the 
effect on the compressibility sum rules has to be exam- 
ined. In any case the DFT implementations in codes like 
VASP, or SIESTA cannot deal correctly with the case 
Ti ^ Te , and it is not clear if they treat the Hi^i term in 
the partition function correctly even in the equilibrium 
case, due to the use of the Born-Oppenheimer approxima- 
tion which prevents the possibility of coupled electron-ion 
plasma modes in the system. 

The CHNC technique is a non-dynamical method that 
does not need the Born-Oppcnhcimcr approximation. It 
correctly treats the cross-interaction H^i even for two- 
temperature systems, as established by direct MD sim- 
ulations 0. Fig. [5] shows that the SESEME and other 
standard EOS agree with the CHNC-BO calculation 
where Tei is set to (Tg -I- Ti)/2, while the Laser-shock 
experiments, where Ti ^ Tg may hold, should agree 
with Tei chosen as the temperature of the scattering pair. 
Ion masses are much larger than mg, and hence T^i ap- 
proaches the electron temperature, as demonstrated in 
Dharma-wardana and Murillo via MD simulations [25| . 
In effect, the calculation of the Laser-shock hydrogen 
Hugoniot has to address non-equilibrium effects, as well 
as non-adiabatic effects associated with the use of the 
Born-Oppenheimer approximation in standard simula- 
tions. The conclusions of Galli et al. [i^l also point to 
non-equilibrium effects associated with the electron-ion 
interaction, i.e., precisely the term iJei in the Hamil- 
tonian indicated in Eq. 1221 Our own views have 
evolved beyond what we stated in Ref. [26|, and the 
subject probably needs to be revisited, within a two- 
temperature quasi-equilibrium setting, without making 
the Born-Oppenheimer approximation, especially at very 
high compressions. 




0.6 O.S 

density (g/cc) 

FIG. 5: Comparison of the CHNC Hugoniot with experiment 
and other theories for warm-dense Deuterium. Of these, the 
SESEME (SES), path-integral Monte-Carlo (PIMC), CHNC- 
BO and the CHNCl are of interest. A non-equilibrium Hugo- 
niot, marked NEQ is also shown [2^. Experiments 1, 2 and 
3 refer to Da Silva et al., Collins et al., and Knudson et al., 
respectively, as described in [2^ . 



Pseudopotentials 

We may also consider the case when the ions are not 
fully ionized into bare nuclei, but carry a group of core 
electrons. For instance, Al-plasmas at 0.5 eV and nor- 
mal compression have a charge Z — i and a core of 10 
electrons. Although it is sufficient for many problems to 
treat the Al"^+ as point charges, a more accurate theory 
may wish to include the effect of the the core radius and 
well-depth of the electron-ion interaction via a pseudopo- 
tential. Such pseudopotentials are well known at zero 
temperature. A very simple model is that of Ashcroft, 
while modern implementations are very sophisticated. 

Al-pseudopotentials suitable for WDM have been given 
in p arametrized from by Perrot and Dharma-wardana 
[49{ . The basic idea is to generate the charge density n(r) 
around a given nucleus of charge Zj^, immersed in a UEF 
of density parameter r^., at a temperature T. The ion is 
place in a spherical cavity in the positive background 
(for details see Ref. [i^) and nir) is determined by a 
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Kohn-Sham calculation which satisfies the Friedel sum 
rule and other properties. Then we define a weak non- 
local pseudopotential Vps {r) by the following relations in 
q-space. 



TABLE I: simple pseudopotentials for Al, C, Si at normal 
compression, and suitable for the WDM regime, in a.u. 



Vpsiq,rs,T, Zm) 

xiq) 



An{q,r,,T,ZN)/x{q,rs,T) (23) 
x{q)V{^-Ml-G,)x{q,r,,T)" 



Here x''(g,rs,r) is the Lindhard response function at fi- 
nite T and electron Wigner-Seitz radius rg, Vq = An/q'^ 
and Gq is a local field correction consistent with the 
density and temperature of the UEF. Further more, 
An(g, rs,T, Zpf) is the Fourier transform of the real-space 
free-electron-density pileup An{r) calculated at the jel- 
lium density n — 3/{ATrr^) and temperature T, for the 
nucleus Z„. That is 



An{r) = nf{r)—n 
nf{r) — n{r) ~ rihir) 



(24) 
(25) 



The bound electron density nb{r) is obtained from the 
orbitals of the finite-T Kohn-Sham equation as in Pcrrot 
[43| . Here it should be noted that the bound electrons 
have to be assigned to a nucleus keeping in mind that 
some bound states are those of 'hopping electrons' which 
form a band of localized states near the continuum [3] ■ 
Equation l23l defines the pseudopotential to be capable of 
recovering the charge-pile up via linear response. Hence 
it has to render a weak potential. It is not very satisfac- 
tory if the resulting pseudopotential proves to be strong. 
However, the method seems to work in most cases. The 
pseudopotential can usually be parametrized (as in an 
Ashcroft empty-core potential), with a core depth Aq and 
a core radius Tc such that 



Vps{r) = -^0, r <rc 
= Z/r, r > Tc 



(26) 
(27) 



This is evidently a very simple form, compared to mod- 
ern, hard, non-local pseudopotentials used in solid-state 
calculations at T = 0. Such modern potentials remove 
the core, but a Kohn-Sham equation has to be solved as 
they are not weak, and cannot be treated using linear re- 
sponse. However, we have found that simple potentials 
as in Eqs. [23ll26l are adequate for even the liquid-metal 
regime close to the melting point, even for metals which 
require non-local pseudopotentials at T = 0. Excellent 
accuracy is obtained if the response functions x{q) ^-re 
calculated for electrons with an effective mass m* spec- 
ified for each case. It is particularly important to note 
that the 'mean ionization', i.e., Z is a parameter which 
appears in the pseudopotential. The Z is also the La- 
grange parameter defining the total charge neutrality of 
the plasma, as discussed in Refs. 0], [S]- A few ex- 
amples of this type of simple pseudopotentials are given 
in Table HI The C and Si pseudopotentials were used 
to generate PDFs of these ionic liquids and compared 



element 


Rws 


Z 


Ao 




m* 


Al 


3.141 


3.0 


0.3701 


0.3054 


0.998 


C 


2.718 


4.0 


0.0 


0.3955 


1.658 


Si 


3.073 


4.0 


0.0 


0.9475 


0.98 



with Car-Parinello simulations in Ref. [46|. Thus these 
pseudo-potentials can be used in the CHNC equations, 
or in the CMMD simulations, to take account of the exis- 
tence of a finite-sized core. Such methods can be used to 
discuss properties of warm dense matter, thus providing 
a complementary approach to the simulations based on 
statistical potentials discussed by Graziani et al. in the 
context of the Cimarron project for simulations of warm 
dense matter [45 1. 



TWO-TEMPERATURE QUASI-EQUILIBRIA 
AND NON-EQUILIBRIUM SYSTEMS. 

When energy is deposited rapidly in matter using laser 
radiation, the electrons absorb the energy directly and 
equilibriate among themselves, achieving a very high elec- 
tron temperature Te. The ion subsystem, at temperature 
Ti, takes much longer to heat up due to the slow temper- 
ature relaxation via the electron-ion interaction. Hence, 
in laser- heated systems, it is common to find Tg > Ti. 
The inverse situation prevails in shock-heated materials 
since the energy of the shock wave couples to the heavy 
ions and not to the electrons [2|. 

The possibility of using a static approach like the 
CHNC for non-equilibrium systems resides on Bogoli- 
ubov's idea of timescales and conserved quantities. We 
have exploited these ideas in our work on hot-electron 
relaxation, both within Green's-function methods, and 
within CHNC methods [49|. The parameters Te, Ti 
in a two-temperature system are merely Lagrange pa- 
rameters which assert that, for certain time scales Tg, 
Ti, the subsystem Hamiltonians H,,, Hi are conserved 
quantities. Similarly, a number of other parameters, 
e.g., quasi-equilibrium chemical potentials, thermody- 
namic potentials, pseudopotentials, Z, etc., attached to 
the subsystems may be conserved for the selected time 
scales. In fact, the original discussions of quasi-equilibria 
by Bogoluibov were used in Zubarev's theory of non- 
equilibrium Green's functions, and RPA-like results for 
the quasi-thermodynamics as well as energy relaxation 
were addressed there-in. However, RPA-like theories 
are of limited value. In strongly coupled regimes, the 
PDFs associated with the given subsystems can be con- 
structed using CHNC, where the use of the correct inter- 
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subsystem temperatures (e.g., Tei) for evaluating inter- 
system PDFS (e.g., gei) is essential. The nature of this 
inter-system temperature is revealed by its appearance 
in the inter-subsystem energy- relaxation formula [i?! . A 
calculation of the distribution functions of two tempera- 
ture plasmas using HNC methods as well as MD methods 
was given recently p5| . 



INHOMOGENEOUS SYSTEMS 

The classical-map technique uses a classical fluid at 
a finite temperature Tcf to represent a uniform-density 
quantum fluid at T = 0. The parameter Tcf is density 
dependent, and hence the extension to a system with 
an inhomogeneous density is not straight-forward. Fur- 
thermore, integral-equation techniques like the HNC be- 
come very complicated when applied to inhomogeneous 
systems. Molecular dynamics can be applied if a viable 
mapping can be constructed. However in this connection 
we should note that studies of confined classical electrons 
in parabolic traps have also yielded useful insights [isj. 

The classical-map technique treats the zeroth-order 
Hamiltonian exactly, i.e., the map is constructed to re- 
produce the known quantum solution classically, requir- 
ing the confining potential to be mapped as well. Even 
when there is no confining potential (other than a uni- 
form background), the zeroth-order problem of has 
to be correctly treated. This was done in the UEF prob- 
lem by constructing a Pauli exclusion potential (i.e., the 
Fermi-hole potential) to recover the Fermi hole in the 
9sAr) exactly. 

When non-interacting electrons are placed in an ex- 
ternal potential, e.g., a parabolic trap, the uniform den- 
sity n modifies to a new distribution rPir). Classically, 
this distribution is of the Boltzmann form, Eq. [2] where 
VKs{r) contains all the terms found in the exponent of 
the HNC equation. Thus, given the vP{r) calculated from 
a quantum mechanical treatment of iJo which contains 
the parabolic external potential, it is necessary to invert 
the HNC equation to get the effective classical poten- 
tial which corresponds to n°(r)._A simplified approach 
to this was used by us in ref. 



24|. At this stage the 



calculation is somewhat similar to the determination of 
the Pauli exclusion potential, and hence the specification 
of an effective fiuid temperature does not become nec- 
essary. The classical Coulomb fiuid at a finite tempera- 
ture Tcf is still necessary for dealing with the many-body 
effects generated from the Coulomb interaction. How- 
ever, given a non-uniform distribution, there is no evi- 
dent method of defining a unique Tcf and the simplic- 
ity of the original CHNC method is lost. Further more, 
the electron-electron pair-distribution functions now de- 
pend explicitly on two coordinates, viz., g{ri,r2)ss' ■ The 
use of molecular-dynamics simulations is more convenient 
in dealing with systems where the simplicity of homoge- 



neous systems is lost. Another advantage an MD simu- 
lation is that the the bridge-function approximations are 
avoided. 

In mapping an inhomogeneous system of density n(r) 
to a hornogeneous slab of density Ti we have used the form 

lilillli. 



:< n{r)n{r) >< n{r) > 



(28) 



in dealing with 2D systems. The same method has been 
used by Gori-Giorgi and Savin for 3D systems 51]. Using 
such a uniform density to define a unique temperature 
of an equivalent classical fluid, we were able to repro- 
duce the charge distribution of interacting electrons in 
2D quantum dots obtained from Quantum-Monte Carlo 
methods. However, as we used CHNC, it was necessary 
to introduce bridge-functions and boundary corrections 
which impaired the transparency of the classical map. 
Hence this work 2M may be regarded as a preliminary 
attempt. 



CONCLUSION 

We have outlined the classical-map technique of treat- 
ing the quantum many-body problem in Fermi systems 
via a mapping to an equivalent classical system at a 
density-dependent effective temperature different from 
the physical temperature, and where the particles inter- 
act by a pair potential consisting of a Pauli-exclusion 
potential and a diffraction-corrected Coulomb potential. 
Large numbers of particles, and their thermodynamics 
or quasi-thermodynamics can be easily calculated. Since 
pair-distribution functions can be calculated accurately, 
and at any value of the coupling constant, the adiabatic 
connection formula provides results for the non-local 
exchange-correlation functionals in an entirely unam- 
biguous, rigorous manner. No gradient corrections, meta- 
functionals etc., are needed. The Born-Oppenheimer ap- 
proximation is not necessary as the CHNC technique is 
not dynamical. Hence the method would be of great in- 
terest from the point of view of equations-of-state studies, 
both for equilibrium, and for quasi-equilibrium systems. 

Since suitable derivatives of the free energy with re- 
spect to density, temperature, and spin polarization lead 
to Landau Fermi-liquid parameters, the method is capa- 
ble of easily furnishing alternative results for the effective 
mass m*, Lande-g factor, local- field factors of response 
functions etc., which are difficult to determine by stan- 
dard Greens-function perturbation techniques of quan- 
tum many-body theory. 

The application of the method to inhomogeneous sys- 
tems is still poorly developed. Similarly, the method, 
being a technique for the total energy as a functional of 
the pair density, is similar to DFT in not yielding spectral 
information within its own formal structure. 
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